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Chapter  1 
Introduction 


The  physical  optics  (PO)  scattered  field  from  perfectly  electrically  conducting  structures  is 
limited  in  accuracy  because  the  PO  current  fails  to  closely  approximate  the  exact  current 
near  surface  discontinuities  such  as  edges.  The  accuracy  of  the  PO  scattered  field  can 
be  significantly  increased  by  adding  the  fringe  wave  (FW)  field,  which  takes  into  account 
the  distortion  of  the  current  caused  by  edges.  An  approximation  to  the  FW  field  can  be 
calculated  by  integrating  physical  theory  of  diffraction  equivalent  edge  currents[l],[2]  along 
the  illuminated  part  of  the  edges  of  the  structure.  These  edge  currents  are  determined  from 
an  integration  of  the  FW  current  (the  exact  current  minus  the  PO  current)  along  incremental 
strips  on  the  canonical  wedge  or  half-plane.  Throughout  this  report  only  physical  theory 
of  diffraction  equivalent  edge  currents  will  be  considered,  and  these  will  be  referred  to  as 
EEC’s. 

Closed-form  expressions  for  EEC’s  have  been  derived  for  un-truncated  (infinite)  incre¬ 
mental  wedge  strips  by  Michaeli  [2],  Mitzner  [3],  and  Shore  and  Yaghjian  [4]^;  these  EEC’s 
will  be  called  un-truncated  EEC’s  in  this  report.  For  the  analysis  of  bistatic  radar  scattering 
the  un-truncated  EEC’s  give  an  inadequate  correction  to  the  PO  field  due  to  the  presence 
of  the  Ufimtsev  singularity  [2]  and  the  discontinuities  of  the  calculated  FW  field  across  the 
current  layers  associated  with  the  un-truncated  strips.  The  Ufimtsev  singularity  occurs  when 
the  direction  of  observation  is  the  continuation  of  an  incident  field  grazing  the  face  of  the 
structure. 

The  above-mentioned  problems  associated  with  the  un-truncated  EEC’s  are  eliminated 
by  using  truncated  (finite)  strips.  Closed-form  expressions  for  EEC’s  have  been  derived  for 
truncated  incremental  half-plane  strips  by  Breinbjerg  [5]  and  by  Shore  and  Yaghjian  [6].  Cote 
et  al.  [7]  have  implicitly  derived  EEC’s  for  truncated  incremental  strips  on  a  right-angled 

Received  for  publication  3  Oct  1995 

^The  results  by  Mitzner  [3]  and  Shore  and  Yaghjian  [4]  are  expressed  in  terms  of  incremental  length 
diffraction  coefficients,  which  are  closely  related  to,  and  can  easily  be  put  in  the  form  of,  EEC’s. 
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wedge.  Michaeli  [8]  seems  to  be  the  only  one  who  has  derived  EEC’s  for  general  truncated 
incremental  wedge  strips;  the  EEC’s  based  on  truncated  wedge  strips  will  in  this  report  be 
called  truncated  EEC’s.  They  apply  to  the  analysis  of  bistatic  radar  scattering  from  three- 
dimensional  structures  with  flat  faces.  However,  from  theoretical  considerations,  as  well  as 
numerical  calculations,  it  appears  that  Michaeli ’s  truncated  EEC’s  contain  non-removable 
singularities  which  give  rise  to  numerical  problems  and  thus  hamper  their  application.  The 
singularities  are  caused  by  the  mathematical  procedure  applied  to  obtain  closed-form  expres¬ 
sions,  and  they  occur  for  special  directions  of  incidence  and  observation  and  for  zero  strip 
length. 

In  this  report  new  truncated  EEC’s  are  derived.  These  EEC’s  do  not  have  the  above- 
mentioned  singularity  problems  of  the  previously  reported  expressions  [8],  that  is,  they  are 
well-behaved  for  all  directions  of  incidence  and  observation  and  they  take  a  finite  value 
for  zero  strip  length.  The  new  truncated  EEC’s  are  thus  well-suited  for  application  to  the 
analysis  of  bistatic  radar  scattering  from  three-dimensional  structures  with  flat  faces. 

The  report  is  organized  as  follows.  In  Chapter  2  the  concept  of  truncated  EEC’s  is 
introduced.  In  Chapter  3  it  is  explained  how  Michaeli’s  truncated  EEC’s  are  derived  and  his 
final  expressions  are  presented.  The  derivations  that  lead  to  the  new  truncated  EEC’s  are 
performed  in  Chapter  4,  and  finally.  Chapter  5  presents  numerical  examples  to  illustrate  the 
differences  among  the  fields  calculated  from  the  method  of  moments  (applied  to  the  magnetic- 
and  electric  field  integral  equations),  the  un-truncated  EEC’s,  Michaeli’s  truncated  EEC’s, 
and  the  new  truncated  EEC’s. 
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Chapter  2 

The  Concept  of  Truncated  EEC’s 

The  configuration  under  consideration  is  a  perfectly  conducting  three-dimensional  structure 
with  flat  faces  illuminated  by  a  plane  wave  (see  Figure  2.1).  In  the  far  field  of  the  structure  a 
high-frequency  approximation  to  the  FW  field  (the  total  scattered  field  minus  the  PO  field) 
is  calculated  from  a  line  integral  along  the  illuminated  part  C  of  the  edges  of  the  structure. 
The  truncated  EEC’s  are  represented  by  the  magnetic  current  Mr  and  the  electric  current 
It  so  that  the  electric  FW  field  is  given  by  [1] 

jk  J  [ZIts  X  (i  X  t)  -|-  Mjs  x  t  ^ ^  dl  (2-1) 

where  j  is  the  imaginary  unit  (the  time  factor  exp(ju?t)  is  suppressed),  k  is  the  wave  number, 
and  Z  is  the  intrinsic  impedance  of  the  ambient  medium.  Furthermore,  s  is  the  unit  vector 
in  the  direction  of  the  far-field  observation  point,  s  is  the  distance  between  the  integration 
and  the  observation  points  and  t  is  the  edge  unit  tangent  vector.  The  two  adjoining  faces 
at  each  edge  are  denoted  by  A  and  B.  It  is  convenient  to  introduce  a  local  rectangular  xyz- 
coordinate  system  with  origin  at  the  integration  point  at  the  edge.  This  coordinate  system 
is  associated  with  face  A;  the  unit  vector  x  is  in  the  plane  of  face  A,  y  is  the  outward  normal 
unit  vector  of  this  face  and  z  equals  the  edge  unit  tangent  vector  t  (see  Figure  2.1).  Using 
this  rectangular  system,  the  vector  s  is  expressed  as  s  =  xsinjS cos  </>  +  y  sin ^ sin ^  -f-  ^ cos ^ 
where  ^  and  <l>  are  the  polar  and  azimuthal  angles,  respectively.  The  propagation  vector  lo  of 
the  incident  plane  wave  is  expressed  as  Sq  =  —x  sin  Po  cos  <l>Q  —  y  sin  sin  ^  cos  where 
/3o  and  <j>o  A  are  the  polar  and  azimuthal  angles,  respectively  (see  Figures  2.1  and  2.2). 
Face  B  is  located  in  the  plane  described  by  ^  =  Ntt  where  Ntt  is  the  exterior  wedge  angle. 
Throughout  the  report  it  is  assumed  that  1  <  N  <2. 

The  truncated  EEC’s  axe  determined  by  a  sum  of  two  contributions,  one  from  each  of 
the  faces  A  and  B 

Mt  =  +  (2.2) 
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Figure  2.1:  Three-dimensional  view  of  a  flat  face  of  a  three-dimensional  structure.  The 
truncated  incremental  strip  extends  from  the  leading  edge  to  the  trailing  edge  and  is  directed 
along  the  unit  vector  The  directions  of  incidence  and  observation  are  denoted  by  io  and 
s,  respectively.  Ntt  is  the  exterior  wedge  angle  and  it  is  assumed  that  1  <  iV  <  2. 


Figure  2.2:  Two-dimensional  view  of  the  configuration  shown  in  Figure  2.1  in  the  plane 
z  =  0. 
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Ir  = 


(2.3) 


Henceforth,  the  superscripts  A  and  B  refer  to  the  contributions  from  the  faces  A  and  B, 
respectively.  In  this  report  the  contribution  from  face  A  will  be  derived  in  detail,  and  the 
contribution  from  face  B  is  then  obtained  from  the  result  for  face  A  using  a  substitution 
technique. 

The  contribution  from  face  A  to  the  truncated  EEC’s,  Mj  and  /y ,  is  calculated  analyt¬ 
ically  by  integrating  the  FW  current  on  face  A  of  a  wedge  appropriately  conforming  to  the 
structure  along  a  truncated  incremental  strip  with  the  length  The  strip  extends  from  the 
leading  edge  (the  edge  at  which  the  EEC’s  are  placed)  to  the  trailing  edge  and  is  directed 
along  the  unit  vector  which  is  the  intersection  of  the  Keller  cone  and  the  face  A,  as  shown 
in  Figure  2.1.  However,  the  integration  of  the  FW  current  along  the  truncated  incremental 
strip  cannot  be  performed  excictly  in  closed  form  and  thus,  an  asymptotic  calculation  is 
necessary.  To  this  end,  the  truncated  EEC’s  are  expressed  as  the  difference  between  the 
un-truncated  EEC’s  and  the  correction  EEC’s, 


Mt  =  MuT  ~  Moot  (2*'^) 

It  =  luT  —  IcoT-  (^•^) 


Michaeli  [2], [8]  found  that  the  un-truncated  EEC’s  can  be  expressed  exactly  in  closed  form 
whereas  closed-form  expressions  for  the  correction  EEC’s  can  only  be  obtained  using  an 
asymptotic  technique.  In  this  report  Michaeli’s  un-truncated  EEC’s  are  used  but  a  new 
asymptotic  calculation  of  the  correction  EEC’s  is  performed  because  Michaeli’s  correction 
EEC’s  contain  singularities  that  give  rise  to  numerical  problems. 

The  contribution  from  face  A  to  the  un-truncated  EEC’s,  and  is  obtained  by 
integrating  the  FW  current  on  face  A  along  an  un-truncated  incremental  strip.  The  strip 
extends  from  the  leading  edge  and  is  directed  along  u^.  Michaeli  found  that  [8,  eqs.  (3)-(7)] 


=  -ZsinPo'^K^ 

smp 

/yy  =  sin  —  sin  0(t  cot  fi  cos  <i>^x 


(2.6) 

(2.7) 


with 

oo 

I^x,z  =  J  exp{jkus  •  du  (2.8) 

0 

where  denotes  the  x-  and  ^-components  of  the  FW  current  on  face  A. 

Similarly,  the  contribution  from  face  A  to  the  correction  EEC’s,  and  is  obtained 
by  integrating  the  FW  current  on  face  A  along  another  un-truncated  incremental  strip.  This 
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strip  extends  from  the  point  of  truncation  at  the  trailing  edge  and  is  directed  along  u^. 
Thus, 


Mi,  = 

Sin  p 

(2.9) 

=  sin  —  sin  cot  ^  cos  <j>L'^ 

(2.10) 

with 

OO 

it  =  / 

(2.11) 

Michaeli  [8]  calls  the  correction  EEC’s  the  second-order  EEC’s  but  has  no  designation  for 
the  truncated  EEC’s.  In  contrast  to  Michaeli’s  approach,  the  correction  EEC’s  in  this  report 
are  placed  at  the  leading  edge  instead  of  the  trailing  edge. 

The  truncated  EEC’s  take  into  account  the  first-order  and  part  of  the  second-order  edge 
diffraction.  The  entire  second-order  diffraction  is  not  taken  into  account  because  the  FW 
current  excited  at  the  trailing  edge  is  neglected.  It  is  possible  to  derive  EEC’s  that  taJke  into 
account  the  FW  current  excited  at  the  trailing  edge  by  employing  the  procedure  introduced 
by  Breinbjerg  [5].  However,  this  is  not  the  concern  of  this  report. 
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Chapter  3 

Previous  Results 


In  this  chapter  the  expressions  for  Michaeli’s  truncated  EEC’s  are  discussed.  To  understand 
the  approach  used  for  the  derivation  of  the  EEC’s,  the  expressions  for  the  exact  FW  current 
on  face  A  of  the  wedge  are  dealt  with  first. 


3.1  The  Exact  Fringe  Wave  Current  on  Face  A  of  the 
Wedge 


The  X-  and  ^-components  of  the  FW  current  on  face  A  of  the  wedge  can  be  expressed  in 
terms  of  a  contour  integral  in  the  complex  ^-plane  [8,  eq.  (9)]  as 


j}w,A 

jfxu,A 


-Hzq  exp{-jkz  cos  ^o)Ix 


cot 

ttN  V  Z  sm  pq  ) 


with 


h 

h 


exp{jX  cos^)d( 
expijX  cos()d( 
exp{jX  cos^)  d^. 


(3.1) 

(3.2) 

(3.3) 

(3.4) 

(3.5) 


In  these  expressions  F  is  the  steepest  descent  path  through  tt  (see  Figure  3.1)  and  X  = 
kxsin^Q.  Furthermore,  Hzo  and  Ezo  are  the  z-components  of  the  incident  magnetic  and 
electric  field,  respectively,  at  the  origin  of  the  local  xyz-coordinate  system. 
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Figure  3.1:  The  steepest  descent  path  F. 

The  exact  FW  current  given  by  (3.1)  and  (3.2)  is  found  in  the  following  way.  First,  the 
total  field  is  found  for  illumination  by  an  electric  and  magnetic  z-directed  line  source.  This 
is  achieved  by  employing  an  expansion  in  terms  of  eigenfunctions  [9,  sec.  3.3].  The  solution 
for  plane  wave  TE  and  TM  (with  respect  to  z)  incidence  is  next  obtained  by  moving  the 
line  source  to  infinity.  The  solution  for  oblique  incidence  is  obtained  by  employing  the 
substitution  technique  described  by  James  [9,  p.l27]  and  Bowman  et  al.  [10].  Next,  the 
solution  for  the  total  field  is  expressed  in  terms  of  contour  integrals  and  the  diffracted  field  is 
finally  arrived  at  by  deforming  the  integration  contour.  Applying  the  result  for  the  diffracted 
field,  the  components  (3.1)  and  (3.2)  of  the  FW  current  are  easily  found. 

3.2  Un- Truncated  EEC’s 

The  calculation  of  the  contribution  from  face  A  to  the  un-truncated  EEC’s  (2.6)  and  (2.7) 
is  accomplished  by  insertion  of  the  expressions  (3.1)  and  (3.2)  for  the  exact  FW  current 
into  the  integral  in  (2.8).  Next,  the  order  of  integration  is  interchanged  and  the  inner 
integral  is  calculated  analytically.  Finally,  the  integration  contour  F  is  distorted  and  closed 
at  infinity  and  thus  the  integral  can  be  calculated  by  applying  the  residue  theorem.  The 
results  are  [2,  eqs.  (4)- (7)], 

_  2jZHzosm<l>  /  U{Tr-(j>o)  _ 

A;  sin  ^0  sin  ^  \ /X  +  cos  <^o  A’sina(cos  —  cos  ^) 
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and 


M  — 
lUT  — 


2j 


JbiVsin^o(cos  ^-coslf)  V  sin  a  v— -  —  Zsin^o 

2jU{'K-(i>o)  ( rr  ,  .  ^  1,  xo  ■£iosin<^o^ 

+  ,  .  ^  ;  .  ,  X  Hzo{coi  pcos<f>  +  cot  ^0  cos  <i)o)  -  y  .  1 

A;  sm  Po{ii  +  cos  cpo)  \  ^  sm  po  / 

2j  cot  ^qHzO 


) 


kN  sin  Po 

where  t/(a;)  =  1  for  rc  >  0  and  zero  for  x  <  0.  Furthermore, 
sin  po  sin  P  cos  (f)  +  cos  ;0o(cos  P  —  cos  Po) 


fj.  = 


sm 


and  a  is  the  solution  to  fj,  —  cos  a  determined  by 


a  =  -jLog{iJ,  +  -  1) 


(3.7) 


(3.8) 


(3.9) 


with 


Logz  =  ln|z|  +  Argz 


(3.10) 


and  —X  <  Argz  <  tt.  The  square  root  in  (3.9)  is  defined  as 


iWi-  ti^\ 


^l  <  —\ 

— 1  <  <  1. 


(3.11) 


The  results  for  and  are  obtained  from  the  above  expressions  by  replacing  Pq  with 
TT  —  Po,  P  with  TT  —  P,  <f>o  with  Ntc  —  <f)o,  and  (f)  with  <f>  —  Ntc.  Using  these  substitutions  it  is 
noted  that  the  last  term  in  (3.7)  is  canceled  by  its  counterpart  in  7®^. 

The  singularities  of  M^j  and  I^x  have  been  carefully  investigated  by  Michaeli  [2].  The 
only  non-removable  singularity  is  the  Ufimtsev  singularity  [2]  which  occurs  when  fi  =  1 
(a  =  0)  and  simultaneously  <j>o  =  Tt,  that  is,  when  the  direction  of  observation  is  the  con¬ 
tinuation  of  an  incident  field  grazing  the  face  of  the  structure.  The  Ufimtsev  singularity  is 
caused  by  the  lack  of  phase  variation  in  the  integrand  of  in  (2.8)  for  large  values  of  the 
integration  parameter  u  and  is  eliminated  by  using  truncated  strips. 


3.3  Michaeli’s  Correction  EEC’s 

The  approach  used  by  Michaeli  [8]  to  calculate  the  contribution  from  face  A  to  the  correction 
EEC’s,  Mt  in  (2.9)  and  in  (2.10),  is  as  follows.  First,  the  expressions  for  the  exact 
FW  current  (3.1)  and  (3.2)  are  inserted  into  the  integral  in  (2.11).  Next,  the  order 
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of  integration  is  interchanged,  the  inner  integral  is  calculated  analytically  and  finally,  the 
resulting  integral  is  evaluated  asymptotically.  However,  there  are  two  problems  associated 
with  this  approach.  First,  the  integrand  of  the  resulting  integral  contains  many  poles  and 
this  fact  complicates  the  task  of  obtaining  an  asymptotic  expression  without  non-removable 
singularities.  In  fact,  Michaeli’s  asymptotic  evaluation  results  in  expressions  that  contain 
non-removable  singularities  for  N  ^2.  These  singularities  occur  for  special  directions  of 
incidence  and  observation.  Second,  the  asymptotic  evaluation  implies  that  when  N  ^2 
the  resulting  expression  has  a  non-removable  singularity  for  L  =  0,  where  L  =  kl^  sin^  /So- 
The  above-mentioned  singularities  give  rise  to  numerical  problems  and  thus  hamper  the 
application  of  the  expressions. 

To  illustrate  the  above-mentioned  problems,  Michaeli’s  asymptotic  evaluation  of  is 
presented  below.  To  distinguish  Michaeli’s  result  from  the  expressions  that  will  be  derived 
in  this  report,  Michaeli’s  result  is  denoted  by 

After  the  interchange  of  the  order  of  integration  and  the  analytical  calculation  of  the 
inner  integral,  the  expression 


» _  ZHzo  sin  <j> 

jNTrksmPo  sin  j3 

(3.12) 

is  arrived  at  where  V  is  the  contour  integral 

Y  ^  f  -j  N  +  cos  0) 

Jr  (cos  jj  -  cos  ^){fi  -1-  cos 

(3.13) 

with  the  contour  F  shown  in  Figure  3.1.  Next,  the  substitution 

/2  ( 

(3.14) 

is  employed  to  transform  the  steepest  descent  path  F  to  the  real 
V  in  (3.13)  becomes 

5-axis  whereby  the  integral 

nr  oo 

v  =  yjexp(jl(/i  -  1))  J  kF(^)exp(-Ls^)d5 

(3.15) 

where  y/J  means  exp(j^)  and 

sin  |(cos  -  cos  ^)(/z  +  cos 

(3.16) 

To  isolate  the  singularities  of  VF(^)  close  to  the  saddle  point  ^  = 
a  decomposition  of  the  form 

TT  (s  =  0),  Michanli  applied 

(3.17) 
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(3.18) 


with  the  assumption  that  5'(^)  is  a  regular  function, 


Ds  +  E 

+jb^' 


In  the  expression  (3.17)  for  ^  ~  y/^cos^  and  b  —  V^sin^.  The  first  term  on 

the  right  side  of  (3.17)  takes  into  account  the  pole  of  (cos  ^  —  cos  ^)~^  near  the  saddle 
point  (^  =  tt)  when  <^o  =  tt.  The  second  term  on  the  right  side  of  (3.17)  takes  into  account 
the  two  poles  of  {fi  +  cos{)“^  near  the  saddle  point  when  a  =  0  (recall  that  n  =  cos  a). 
There  is  a  problem  associated  with  the  decomposition  (3.17).  The  problem  occurs  when  the 
decomposition  constant  E  is  to  be  determined^.  To  calculate  this  constant  one  multiplies 
equation  (3.17)  by  +  jb"^  on  both  sides.  It  is  seen  that  fi  +  cos  ^  >  0  and  +  jb"^  0  for 
^  IT  ±  a.  Thus,  by  employing  the  fact  that  5'(^)  is  regular  and  by  calculating  the  limits 
for  ^  >  tt  —  a  and  ^  »  tt  +  a,  and  combining  the  two  results,  the  constant  E  is  determined. 

However,  if  JV  7^  2  and  a  =  ir  the  quantity  sin  ^(sin  |)”^  appearing  in  (3.16)  tends  to  infinity 
for  ^  — V  z  +  a  and  this  singularity  is  not  accounted  for  by  the  decomposition  (3.17).  Since 
S'(^)  is  assumed  regular  the  singularity  for  a  =  tt  will  appear  in  the  decomposition  constant 
E.  For  iV  =  2  the  quantity  sin^(sin|)“^  equals  1,  and  no  problem  occurs.  Further,  for 
(f)Q  =  -(ff  +  a)  +  2TrN  the  quantity  (cos  jf  —  cos  ^)“^,  which  appears  in  the  expression  (3.16) 
for  W(^),  goes  to  infinity  for  ^  x  +  a.  If  iV  =  2  this  singularity  will  be  accounted  for  by 
the  term  C{s  +  ^)"^  in  (3.17).  If,  on  the  other  hand,  N  ^2  there  will  be  no  term  on 
the  right  side  of  (3.17)  that  will  take  this  singularity  into  account  and  thus,  the  singularity 
appears  in  E.  In  summary,  the  decomposition  constant  E  contains  two  singularities  when 
N  ^2.  These  occur  for  a  =  tt  and  for  ^0  =  “(’t  +  «)  +  27riV. 

The  first  two  terms  on  the  right  side  of  (3.17)  are  integrated  explicitly  using  (3.15)  and 
expressed  in  terms  of  Fresnel  functions.  The  third  term  of  (3.17)  is  integrated  asymptotically 
for  £  >  1  using  the  standard  steepest  descent  technique  [11].  For  N  ^2  this  technique 
implies  that  the  resulting  expression  tends  to  infinity  for  L  0.  This  can  occur  for  edge 
points  close  to  corners  in  the  evaluation  of  the  integral  (2.1).  As  discussed  by  Michaeli  [8]  it 
is  possible  to  avoid  small  values  of  L  by  omitting  part  of  the  edge  which  is  close  to  corners 
in  the  evaluation  of  the  integral  (2.1).  However,  this  approach  is  not  satisfactory  since  the 
calculated  field  will  depend  on  the  ratio  of  the  edge  being  omitted.  For  N  —2  the  asymptotic 
calculation  equals  the  exact  solution  and  thus,  no  singularity  occurs  when  £  =  0. 

To  summarize,  the  singularities  occurring  for  <j)Q  =  —  (tt  +  q;)  +27rA^  and  a  =  tt  are  caused 
by  the  decomposition  of  the  integrand  (3.16)  that  contains  many  poles.  It  is  noted  that 
the  quantity  fi  +  cos^  in  the  denominator  of  (3.16)  is  caused  by  the  analytical  calculation 
of  the  inner  integral  after  the  order  of  integration  is  interchanged  (see  the  summary  of 

^It  is  not  necessary  to  determine  the  constant  D  because  the  odd  term  vanishes  when  integrated 

from  —00  to  00,  see  (3.15). 
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Michaeli’s  approach  at  the  beginning  of  this  section).  Thus,  the  interchange  of  the  order 
of  integration  and  the  following  analytical  calculation  of  the  inner  integral  introduce  more 
poles  to  the  integrand.  Furthermore,  the  singularity  for  T  =  0  is  caused  by  the  use  of  the 
standard  steepest  descent  technique  on  the  integral  (3.15)  containing  the  exponential  factor 
exp(— Z/5^).  This  exponential  factor  is  introduced  by  the  analytical  calculation  of  the  inner 
integral.  It  is  therefore  seen  that  the  singularity  problems  of  Michaeli’s  expressions  are 
basically  caused  by  the  interchange  of  the  order  of  integration  when  the  integral  in 
(2.11)  is  calculated.  In  the  following  chapter  the  integral  will  be  calculated  without 
interchanging  the  order  of  integration.  By  this  approach  the  singularity  problems  described 
in  this  section  are  avoided. 

The  asymptotic  results  for  M^j.  and  valid  for  >  1,  are  [8] 


= 


ZHzo  sin  <f>exp{jL{fi  -  1)) 
j%N  k  sin  ^0  sin  /3 


fJL  +  COS  00  ^ 


^/__sin^_^ 
sin  a  \cos  —  cos  ^  ( 

1^ _ sin^ _ 

V  jL  (cos  ^  -  cos  ^)(fi  -  1)J 


N  N  / 


(3.19) 


= 


_  2exj>{jL{fi  -  1))  f  A'^sgn(cos  ^) 
jkN  sin  \  ^  +  cos 


•  ( -ii4o(cot  ficos<t>  +  cot  /3o  cos  ^o)  +  K{V^\  cos  ^|) 

^  sin  fjQ  J  Z 


+ 


'  (  sin^ 

sin  a  \cos  2^  —  cc 

sin  ^  / _ 1 

Zsin;0o  Vcos^^- 


<f>0 

N 


cos^-cos^/ 


cos 


<i>o 

N 


cos  ^  -  cos  ^ 


{fi  cot  |5o  —  cot  ^  cos  <l>) 


Hzo  sin  ^(cot  ^  cos  <f>  —  cot 


y/2'KjL{cOS  jj-COSjf 
In  these  expressions  sgn(a)  =  ^  and 

K{^)  =  F{^)  -  2^^ 

where  F  is  a  modified  Fresnel  function 
F{x)  =  \J^exi[>{jx^)  J  exp(-jt^)dt. 


I 

-l)i‘ 


(3.20) 


(3.21) 


(3.22) 
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It  is  noted  that  the  expressions  for  and  contain  non-removable  singularities  when 
N  ^2.  These  singularities  occur  for  T  =  0,  a  =  tt,  and  <j>Q  =  — (tt  -|-  a)  -|-  2TtN. 

The  expressions  for  and  7®,  are  obtained  by  applying  the  substitutions  listed  after 
equation  (3.11).  Furthermore,  the  strip  length  on  face  A,  must  be  replaced  by  the  strip 
length  on  face  B,  l^. 
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Chapter  4 

Derivation  of  New  Correction  EEC’s 


As  summarized  in  Section  3.3  the  singularity  problems  associated  with  Michaeli’s  correction 
EEC’s  are  caused  by  the  interchange  of  the  order  of  integration  when  calculating  the  integral 
Li  g  in  (2.11).  In  this  chapter  another  approach  is  used  to  calculate  Instead  of  applying 
the  exact  expressions  for  the  FW  current  in  the  integral  the  asymptotic  expressions 
for  the  FW  current  are  employed.  Thus,  the  first  section  in  this  chapter  deals  with  the 
determination  of  the  asymptotic  expressions  for  the  FW  current. 


4.1  Uniform  Asymptotic  Expressions  for  the  Fringe 
Wave  Current  on  Face  A 


To  obtain  a  uniform  asymptotic  expression  for  the  FW  current,  the  steepest  descent  tech¬ 
nique  [11]  is  applied  to  the  integrals  hyh,  and  I3  given  by  (3.3)-(3.5).  The  steepest  descent 
path  r  is  first  transformed  to  the  real  s-axis  using  the  substitution  (3.14) 


where  y/J  means  exp(j|).  For  Ii  this  yields 


00 

h  =  y/^exp{-jX)  j  T{i)ex^{-Xs‘^)ds 

—00 


(4.1) 


(4.2) 


where 


m  = 


Sin 


N 


(4.3) 
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To  isolate  the  pole  of  the  integrand  near  the  saddle  point  (^  =  x)  a  decomposition  technique 
similar  to  the  one  applied  by  Michaeli  [8,  eq.  (29)]  is  used.  Thus,  T{()  in  (4.3)  is  written  as 
a  sum  of  a  simple  pole  term  and  a  regular  term, 


rK)  = 


s+  — 


+m) 


(4.4) 


where  a  =  ^/2  cos  ^  and  the  regular  term  is 


m) = no  - 


(4.5) 


The  poles  of  the  first  term  on  the  right  side  of  (4.4)  are  determined  by  ^  =  i^o  +  47rn,  n 
being  an  integer,  and  thus,  the  pole  at  ^  =  x  is  isolated  for  <j>o  =  x.  The  decomposition  (4.4) 
contains  only  one  simple  pole  term  and  is  thus  simpler  than  the  decomposition  presented 
in  Section  3.3.  Due  to  the  simpler  decomposition,  the  singularity  problems  encountered  in 
Section  3.3  will  not  appear. 

The  decomposition  constant  A  is  determined  from 


yf  sin^(-cosf +  cos^)  _  ^ 
sin  |(cos  jf  —  cos  ^) 

The  integration  of  the  first  term  on  the  right  side  of  (4.4)  is  written  as 


(4.6) 


7  exp(-££)  ^^7 


(4.7) 


The  integral  on  the  right  side  of  (4.7)  is  then  evaluated  using  [12,  eq.  (12)]  and  expressed  in 
terms  of  the  modified  Fresnel  function  (3.22) 


F{x)  =  ^exp{jxO  J  exp{-jtOdt 


(4.8) 


with  the  result 


/  ds  =  2ixsgn(a)F(|a|VX) 


(4.9) 


where  sgn(a)  =  Since  the  function  R{^)  is  regular  near  the  saddle  point,  the  standard 
steepest  descent  technique  is  applied  to  the  second  term  on  the  right  side  of  (4.4), 

oo  - 

J  R{Oexp{-XsOds^  R{x)>Jj  for  X  >  1.  (4.10) 

—  OO 
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Thus,  the  asymptotic  expression  for  Ii  is 

h  ~  -2j7r7Vexp(-jA')sgn(a)F(|<i|A/^) 


1^  /  sin^  N\ 

V  IT  +  'IT )  • 

V  \cos  —  cos  ^  V2a / 


(4.11) 


The  same  procedure  is  used  to  calculate  the  asymptotic  expressions  for  I2  and  I3.  The  results 
are  (see  Appendix  A) 


—2jNir  sin  <^o  exp(— j  A") 


sin 


sgn( 


N 


a)  (fMVX) 


1 


2y/jTrX\a\ 


(4.12) 


—2jirNsgn(a)cos  ^oexp(— jX)F(|a|\/X) 

V  ^  V  V  2a  cos  ^  —  cos  ^  ^ 


(4.13) 


The  asymptotic  expressions  for  in  (3.1)  and  in  (3.2),  valid  for  X  =  fcxsin  ^0  >  1, 
are  obtained  by  employing  the  above  results  for  /i,  I2,  and  h, 


jfva,A 


—2Hzoexp{—jku‘^  •  r) 

/  f  > 

y/2j'Kkx  sin  \2co£ 


sgn(cos  ^)F{^2kx  sin/3o|  cos  ^|) 


+ 


•  TT 

si^ 


cos  ^  iV(cos  ^  —  cos  ^ 


(4.14) 


and 


jSw,A  ^  _2  exp(-j ku^  •  r)  ( -  H^o  cot  ^0  cos  <f)o 

L  \  Z  sm  po  } 

•sgn(cos  ^)F{‘^2kx  sin  Po\  cos  — |) 


Ezo  sin  <f>o 


Z  sin  /9o2  v^JttTxsuT^  cos  ^ 
HzQ  cot  ^0  (  cos  ^0 


(gp  / 

sin^o  V 


sin^ 


V'2j7rA:xsin/3o  \2  cos  ^  N{cos  jj-cos^ 


(4.15) 


where  f  =  xx  +  and  «^  =  x  sin  l3o  +  z  cos  Pq.  The  formulas  for  the  FW  current  can  also 
be  expressed  in  terms  of  the  function  K{x)  in  (3.21), 


K{x)  =  F{x) 


1 


2y/pKX 


(4.16) 
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where  K(x)  «  0  for  x  >  Z.  One  finds 


jfw,A 


-2Hzo  ex^{—jku^  •  r  )  |^sgn(cos  ^)K{\j2kxsm  f3o\ cos  ^  j) 

\ 

3m^(cos  ^  -  cos  j^)) 


Ny/2jxkxsin  i 


(4.17) 


— 2exp(— •  f 


sin^o 

siny^o 


—  Hzo  cot  ^0  cos  I 


•) 


•sgn(cos  ^)K(‘\j2kx  siny^ol  cos  ^|) 


Hzo  cot  /?o  sin  jj 


N y/2jincx^n^{cos  ^  —  cos 


(4.18) 


The  representations  (4.14)  and  (4.15)  are  well-suited  for  the  analytical  calculation  of  the  FW 
current  in  the  limit  4>o  tt,  whereas  using  the  equations  (4.17)  and  (4.18)  is  more  convenient 
to  obtain  expressions  for  the  FW  current  when  ‘\/2kx sin^ol  cos  ^|  >  3. 


4.2  Calculation  of  the  Integrals  and  Lf 


To  derive  the  results  for  M^j.  in  (2.9)  and  7^,  in  (2.10),  the  expressions  for  7/^  and  given 
in  (2.11)  are  needed.  First,  is  calculated.  To  this  end  the  ^-component  of  the  asymptotic 
FW  current  (4.14)  is  inserted  into  the  integral  L^.  The  integration  variable  u  of  is  the 
distance  along  the  Keller  cone  and  thus  the  substitutions  a:  =  u  sin  0o  and  z  —  u  cos  0q  axe 
applied.  Since  the  asymptotic  expression  for  the  FW  current  applies  for  ku  sin^  0q'^  \  the 
asymptotic  result  for  applies  for  L  —  kl^  sin^  0q'^ 


L 


A 

X 


+ 


oo 

J  exp{jkus  -u^)du 

— 2.fl^osgn(cos  -^)  J  F{\/2ku  sin  0q\  cos  -^|)  exp{—jku  sin^  0o{\  —  fi))  du 

27^0  /  1  sin^  \  Jexp{-jku  sin^  0o(l  -  fi)) 

v/2j¥Fsin/3o  \2cos  ^  W(cos  ^  —  cos  ^)  y  ^7  >/u 

(4.19) 
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where  the  identity  1  —  s  •  =  sin^  ^o{\  —  fi),  with  fi  defined  in  (3.8),  has  been  applied.  The 

two  integrals  in  (4.19)  axe  evaluated  in  Appendix  B.  Hence, 


2exp(jL(jn  -  l))Ha) 
jk  sin^  Po 

\/5|  cos 


sgn(cos  ^) 


^(F(v^|cos:^l) 


s/2 


1 

+ 


/i  +  cos^io  V  ' 

‘V(v/i(l-ri)) 


2 


y/\  —  H  \2cos  “  N{cos 


(4.20) 


Second,  Lf  is  calculated.  The  asymptotic  expression  for  the  ^-component  of  the  FW  current 
(4.15)  is  inserted  into  the  integral  Lf  in  (2.11)  to  yield 

OO 

Lf  =  J  exp(jkus  •  u'^)du 

~  -2  ( -  Hzo  cot  cos  sgn(cos  ^ ) 

\  ZsmpQ  J  2 

•  J  F{V2ku  sin  fio\ cos  -^1)  exp{—jku  sin^  ^o(l  —  /x))  du 


+ 


Y-^o^ 
\  Z  sii 


sin  <j)Q 


-  HzO  cot  /?o  cos  ^0 


sin^o 

2Hzo  cot  Po  sin  ^ 


■) 


y/2j'Kk  sin  /do  cos  ^ 


Ny/2j'Kk  sin  /do(cos  ^  —  cos  ^)\ 

^  exp{—jkusm^  ^o{l  —  p,)) 


7 


y/u 


du. 


When  the  results  of  Appendix  B  are  applied,  Lf  is  expressed  as 


2 exp{jL{n  Ezo  sin (f>o  „  ^  j  \  sgn(cos 

}  \  ty  •  a  -^zO  cot  Po  cos  ^0 

[  \  Z  Sin  Po 


(4.21) 


fi  +  cos  <f>o  2 


+ 


jksiir  /3o 
Hzo  cot  /do 


( 


cos  ^  cos  ^0  cos  (^0 


+ 


Sin  ^ 


fi  +  cos  <f>o  2  cos  ^  N{cos  ^  —  cos  ^ 


Ezosm^{fi  -1)  '\V2F{y/L{l-fi)) 
Z  sin  ^o{fi  +  cos  ^o)  J  \/l  “  /^ 


(4.22) 
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4.3  Expressions  for  the  Correction  EEC’s 


The  asymptotic  expressions  for  the  contribution  from  faee  A  to  the  correction  EEC’s, 
in  (2.9)  and  in  (2.10),  are  now  obtained  by  employing  the  results  for  in  (4.20)  and 
Lj  in  (4.22), 


Mi 


2ZHzo  sin  (l>exp{jL{fi  —  1)) 

-sgn(cos  ^) 

jk  sin  /3  sin  j3o 

fi  +  cos  <f>o 

+ 


( 


a/1 


V2sin^ 


y/2{n  +  cos  ^o)  cos  ^  N-s/1  —  fi(cos  ^  —  cos  ^) 


m)) 


(4.23) 


2exp(jL(/i  -  1)) 
j  k  sin  l3o{ii  +  cos  ^o)  [ 
f  Ezq  sin  <j>o 


sgn(cos  y  )F(V5ri 


cos  y  I) 


Z  sin  pQ 


—  /4o(cot  /3o  cos  (j)o  +  cot  ^  cos  <f)) 


Eio  sin  ^  ,  Hi 


V2(i-/<)n\/i(i -/<))(- 

Hzo  sin  ;^(/i  +  cos  <^o)(cot  Pq  —  cot  ^  cos  (f>) ' 


Z  sin  ySo  2  cos 


iV(cos  ^ 


cos  ^)(1 


/^) 


■^(cot  ydo  cos  <f>0  +  cot  ^  cos  ^) 
2 

(4.24) 


The  final  expressions  for  the  truncated  EEC’s  are  obtained  by  first  calculating  and 
by  subtracting  the  above  results,  (4.23)  and  (4.24),  from  the  un-truncated  EEC’s,  (3.6)  and 
(3.7),  as  shown  in  (2.4)  and  (2.5).  Second,  the  contribution  from  face  B  is  calculated  using 
the  results  for  and  Jy  by  replacing  with  x  —  ^o,  ^  with  x  —  /?,  with  Nx  —  <j)o,  <f> 
with  Nx  —  <j),  and  with  /®.  Third,  the  contributions  from  the  two  faces  are  added  in  order 
to  determine  Mt  and  Ij,  see  (2.2)  and  (2.3).  Finally,  these  expressions  are  inserted  into  the 
radiation  integral  (2.1)  to  determine  the  approximate  FW  field  from  the  truncated  EEC’s. 

It  is  noted  that  in  (4.23)  and  7^,  in  (4.24)  do  not  contain  singularities  for  a  =  x 
{fi  =  —1),  <l>o  =  — (x  +  a)  +  2xN,  and  X  =  0  as  do  the  previously  reported  expressions  (3.19) 
and  (3.20).  Using  the  result  of  [2,  app.  II]  it  is  shown  that  when  <f>o  ^  x  the  quantities 
sin^(l  —  and  (cot;3o  —  cot/?cos^)(l  —  which  appear  in  the  expressions  for  M^j. 
in  (4.23)  and  7^^  in  (4.24),  remain  bounded  for  — >  1.  li  (l>o  =  x  and  fi  1,  M^j.  and  7^,. 
are  singular  but  this  singularity  (the  Ufimtsev  singularity)  is  canceled  by  the  singularity  in 
Myj-  and  7yj  given  by  (3.6)  and  (3.7),  respectively.  This  means  that  Mj  and  It  are  valid 
for  all  directions  of  incidence  and  observation.  Besides,  the  fact  that  and  7^^  are  finite 
for  X  =  0  implies  that  no  numerical  problems  arise  when  the  strip  length  becomes  small. 
This  can  occur  for  edge  points  close  to  corners  of  the  structure.  However,  the  field  calculated 
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from  the  truncated  EEC’s  for  these  edge  points  is  a  poor  approximation  to  the  exaet  field 
for  two  reasons:  First,  the  asymptotic  expressions  for  the  correction  EEC’s  are  invalid  and 
second,  the  distortion  of  the  current  on  the  structure  near  corners  is  not  taken  into  account. 


4.4  Special  Cases 


If  it  is  assumed  that  the  arguments  \/5L|cos^|  and  y^Zr(l  —  fi)  of  the  modified  Fresnel 
functions  in  (4.23)  and  (4.24)  are  so  large  that  the  asymptotic  formula  [8,  eq.  (40)] 


can  be  applied,  M^j.  in  (4.23)  and  in  (4.24)  become 


—y/2ZHzo  sin  <j)  sin  ^  exp{jL{n  —  1)) 
jk sin  sin —  /i)(cos  ^  —  cos 
y/2HzQ  sin  ^(cot  —  cot  ^  cos  ^)  exp(jZ-(/x  —  1)) 
jksmPQy/j‘KLN{l  -  /z)(cos  -  cos  ^) 


(4.25) 


(4.26) 

(4.27) 


Except  for  <^o  =  +  Q:)  +  2TrN  and  a  =  tt  (//  =  —  1)  these  expressions  are  the  same  as 

those  obtained  using  Michaeli’s  correction  EEC’s  [8]. 

For  the  half-plane,  that  is,  N  =  2,  the  expressions  for  M^j.  in  (4.23)  and  7^,  in  (4.24) 
simplify.  Besides,  the  uniform  asymptotic  FW  current  found  in  Section  4.1  equals  the  exact 
FW  current.  This  means  that  in  (4.23)  and  7^^  in  (4.24)  are  valid  for  any  value  of  L. 
The  correction  EEC’s,  that  is,  the  sum  of  the  contributions  from  faces  A  and  B,  become 


Mr^r  = 


^  4Zgo  sin^eMjI(^  -  1)) 

jA:sm;9sm^o(/^ +  cos^o)  V  '  2  ^  ^  2'^ 


y/2cos  ^ 

-v/r=7  ^ 


yui  -  /<))) 


(4.28) 
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and 


(21  cos  ^|F(v^| 

jA:Zsm^;So(// +  cos^o)  V  2 

-^2(1  -  ^)F(^L(\  -  ,,))) 

4g^exp(j£(M-l))/ 
j A:  sin  ;0o(Ai  +  cos  ^o)  V  2 


cos 


•(cot  /?ocos  4>q  +  cot  ^  cos  (f)F{y/^\  cos  ^|) 

Zi 

y/2  cos  ^(cot  ^cos(f>-n  cot  ygp) 


(4.29) 


which  are  the  same  results  obtained  using  Michaeli’s  expressions  [8].  The  expressions  (4.28) 
and  (4.29)  are  further  verified  by  letting  =  0.  In  this  case  Meat  =  Mut  and  /cor  =  Iut- 
Thus,  the  truncated  EEC’s,  Mt  in  (2.4)  and  It  in  (2.5),  are  zero  which  is  the  result  obtained 
by  integrating  the  FW  current  along  a  strip  with  length  zero. 
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Chapter  5 

Numerical  Results 


In  this  chapter  the  field  scattered  by  a  two-dimensional  perfectly  conducting  triangular  cylin¬ 
der  is  calculated  using  both  the  EEC’s  and  the  method  of  moments  applied  to  the  magnetic- 
and  electric  field  integral  equations  (MFIE  and  EFIE,  respectively).  The  consideration  of  a 
two-dimensional  structure  is  convenient  because  it  has  no  corners.  The  purpose  of  perform¬ 
ing  the  numerical  comparison  is  twofold.  First,  it  is  shown  that  the  un-truncated  EEC’s  are 
inadequate  for  bistatic  analysis.  Second,  the  singularity  problems  associated  with  Michaeli’s 
expressions  for  the  truncated  EEC’s  are  illustrated  and  it  is  shown  that  these  singularities 
are  not  present  in  the  new  truncated  EEC’s  derived  in  this  report. 

Both  the  FW  field  and  the  total  scattered  field  are  calculated.  By  considering  the  FW 
field  it  is  possible  to  see  how  well  the  EEC’s  approximate  the  exact  field  in  regions  where 
the  PO  field  is  large  compared  to  the  FW  field.  The  FW  field  is,  however,  not  used  without 
the  PO  field  in  practice.  Therefore,  the  total  scattered  field  is  also  calculated  to  illustrate  a 
practical  application  of  the  EEC’s. 


5.1  The  Fringe  Wave  Field 

The  lengths  of  the  three  sides  of  the  triangular  cylinder  are  all  equal  to  2A,  A  being  the 
wavelength,  and  the  illuminating  field  is  a  transverse  electric  (TE)  polarized  plane  wave 
with  direction  of  incidence  shown  in  Figure  5.1.  The  FW  field  is  calculated  in  the  far  field 
of  the  structure  and  expressed  in  terms  of  the  two-dimensional  radar  cross  section  (RCS). 
The  direction  to  the  far  field  observation  point  is  determined  by  the  angle  (p.  To  avoid  the 
removable  singularities  in  the  EEC’s,  (p  is  assigned  the  values  O.l-t-n  degrees,  n  being  an 
integer  ranging  from  0  to  359.  Figure  5.2  shows  the  FW  field  calculated  from  the  difference 
between  the  MFIE  solution  and  the  PO  solution,  and  from  the  un-truncated  EEC’s.  It  is 
seen  that  the  un-truncated  EEC’s  yield  a  poor  approximation  to  the  exact  scattered  FW 
field:  The  Ufimtsev  singularity  occurs  for  y>=60  degrees  and  the  field  is  discontinuous  across 
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the  current  layers  located  at  </?=120,  180,  240,  and  360  degrees.  In  addition,  discontinuities 
at  (^  =  60  and  ip  =  300  degrees  exist  but  these  cannot  be  seen  on  the  RCS  plot  of  Fig¬ 
ure  5.2.  However,  the  phase  of  the  scattered  field  reveals  the  discontinuities.  Theoretically, 
the  Ufimtsev  singularity  should  only  appear  for  grazing  incidence.  However,  due  to  the  fi¬ 
nite  numerical  accuracy  the  singularity  also  appears  for  nearly  grazing  incidence,  as  in  the 
configuration  under  consideration  (see  Figure  5.1). 

Figures  5.3  and  5.4  show  the  results  when  the  truncated  EEC’s  are  used.  Figure  5.3  shows 
Michaeli’s  results  and  in  Figure  5.4  the  results  obtained  from  the  new  truncated  EEC’s  are 
shown.  From  both  figures  it  is  seen  that  the  Ufimtsev  singularity  and  the  discontinuities 
across  the  current  layers  disappear.  However,  Figure  5.3  reveals  that  five  spikes  occur  in 
the  far  field  obtained  from  Michaeli’s  truncated  EEC’s  when  yj=61,  179,  299,  300,  and  301 
degrees.  These  spikes  are  caused  by  the  non- removable  singularities  in  Michaeli’s  expressions. 
The  spikes  at  9?=61,  179,  299,  and  301  degrees  are  caused  by  the  singularity  occurring  when 
^0  =  —(a  +  tt)  -h  2xA^  (see  the  discussion  in  Section  3.3)  which  is  almost  satisfied  for  edges  B 
and  C  (see  Figure  5.1).  The  spike  at  (p  =300  occurs  because  a  is  close  to  tt  (see  Section  3.3) 
for  edge  B.  As  noticed  from  Figure  5.4  no  spikes  occur  when  the  new  truncated  EEC’s 
are  used.  However,  the  agreement  between  the  two  methods  of  calculation  is  not  perfect. 
This  discrepancy  occurs  because  the  truncated  EEC’s  only  take  into  account  part  of  the 
second-order  edge  diffraction,  as  explained  at  the  end  of  Chapter  2,  and  for  TE  polarization 
the  multiple  interactions  between  the  edges  are  significant. 

Next,  the  incident  plane  wave  is  assumed  to  be  transverse  magnetic  (TM)  polarized  (see 
Figure  5.5).  Figures  5.6,  5.7,  and  5.8  show  the  comparison  among  the  results  obtained  from 
the  difference  between  the  EFIE  solution  and  the  PO  solution,  the  un- truncated  EEC’s, 
and  the  truncated  EEC’s.  For  TM  polarization  the  electric  field  has  only  a  component 
parallel  to  the  current  layers  and  thus  the  field  obtained  from  the  un-truncated  EEC’s  is 
not  discontinuous  across  these  layers.  However,  the  derivative  of  the  field  with  respect  to  <p 
is  discontinuous.  Figure  5.8  shows  that  the  spikes  encountered  when  Michaeli’s  truncated 
EEC’s  are  applied  (see  Figure  5.7)  are  not  present  when  the  new  truncated  EEC’s  are  used. 
Furthermore,  almost  perfect  agreement  between  the  two  methods  of  calculation  is  observed 
in  Figure  5.8.  This  is  because  the  multiple  interactions  between  the  edges  are  weak  for  TM 
polarization. 

In  conclusion,  the  un-truncated  EEC’s  are  inadequate  for  analysis  of  bistatic  scattering 
and  the  numerical  comparison  reveals  that  the  new  truncated  EEC’s,  derived  in  this  report, 
give  a  good  approximation  to  the  exact  scattered  FW  field  even  when  the  distance  between 
the  edges  is  as  small  as  2A. 
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5.2  The  Total  Scattered  Field 

First,  the  configuration  shown  in  Figure  5.1,  where  the  incident  plane  wave  is  TE  polarized, 
is  considered.  Figure  5.9  shows  the  total  scattered  field  calculated  from  the  MFIE  solution 
and  from  the  sum  of  the  solutions  obtained  using  PO  and  the  un-truncated  EEC’s.  It  is 
seen  that  the  un-truncated  EEC’s  yield  a  poor  approximation  to  the  exact  scattered  field. 
In  Figure  5.10  the  sum  of  the  PO  field  and  the  field  obtained  from  the  new  truncated  EEC’s 
is  shown.  The  agreement  with  the  MFIE  solution  is  much  better  than  that  produced  by  the 
un-truncated  EEC’s.  Moreover,  the  agreement  between  the  results  obtained  using  the  new 
truncated  EEC’s  and  the  MFIE  solution  is  much  better  for  the  total  scattered  field  than  for 
the  FW  field  alone  shown  in  Figure  5.4. 

Figures  5.11  and  5.12  show  the  total  scattered  field  when  the  incident  plane  wave  is 
TM  polarized,  as  shown  in  Figure  5.5.  In  Figure  5.11  the  un-truncated  EEC’s  are  used  to 
approximate  the  FW  part  of  the  total  scattered  field  and  again  a  poor  approximation  to  the 
exact  field  is  obtained.  However,  when  the  new  truncated  EEC’s  are  used,  almost  perfect 
agreement  is  achieved  (see  Figure  5.12). 
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RCS/?l  (dB) 


Figure  5.1:  Cross-section  of  triangular  cylinder  with  side  length  2A  illuminated  by  a  TE 
plane  wave. 
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Figure  5.3:  FW  field  for  the  configuration  shown  in  Figure  5.1.  Michaeli’s  expressions  are 
used  to  calculate  the  truncated  EEC’s. 


Figure  5.4:  FW  field  for  the  configuration  shown  in  Figure  5.1.  The  new  truncated  EEC’s 
are  employed. 
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RCS/X  (dB) 


Figure  5.5:  Cross-section  of  triangular  cylinder  with  side  length  2A  illuminated  by  a 
plane  wave. 
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Figure  5.7:  FW  field  for  the  configuration  shown  in  Figure  5.5.  Michaeli’s  expressions  are 
used  to  calculate  the  truncated  EEC’s. 
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Figure  5.8:  FW  field  for  the  configuration  shown  in  Figure  5.5.  The  new  truncated  EEC’s 
are  employed. 
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Figure  5.9:  Total  scattered  field  for  the  configuration  shown  in  Figure  5.1 


Chapter  6 
Conclusions 


New  closed-form  uniform  expressions  for  physical  theory  of  diffraction  equivalent  edge  cur¬ 
rents  have  been  derived  for  truncated  incremental  wedge  strips.  The  new  truncated  EEC’s 
are  well-behaved  for  all  directions  of  incidence  and  observation.  The  expressions  are  asymp¬ 
totic  for  L  1,  L  being  a  parameter  proportional  to  the  strip  length;  however,  they  take  a 
finite  value  when  L  is  zero.  This  implies  that  the  new  truncated  EEC’s  are  well-suited  for  the 
analysis  of  bistatic  radar  scattering  from  perfectly  conducting  three-dimensional  structures 
with  plane  faces. 

Numerical  results  were  presented  showing  that  the  use  of  truncated  strips  is  necessary 
for  the  analysis  of  bistatic  radar  scattering.  Furthermore,  it  was  shown  that  the  singularities 
occuring  in  the  previously  reported  truncated  EEC’s  are  not  present  in  the  new  truncated 
EEC’s  derived  in  this  report. 
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Appendix  A 

Asymptotic  Expressions  for  I2  and  /3 


Using  the  same  procedure  that  was  used  for  the  asymptotic  calculation  of  Ii  in  Section  4.1, 
I2  is  expressed  as 

h  =  exp(-jA)  (^2  J  +  J  i?2(0exp(-Xs2)ds  J  (A.l) 

\  —00  n/j  —00  / 


_  sin^yf(-cos|  +  cos^)  ^  -A  sin  4 
^  sin  |(cos  ^  —  cos  y/^ sin  ^ 

and 

p  /.X  ^  sin^ _ ^ 

(cos  ^ -cos  If)  sm  I  «  +  v7 
The  asymptotic  result  for  A  ^  1  is 

^  -2;A?rsin^oexp(-iA) _ ,_yfipn,/yy  1  \ 

- ^ - =«"('■)  -  2^\a\)  ■ 

The  expression  for  I3  is 


As  =  lim 


cos  ^  sin  cos  |  +  cos 


^0  sin  I  (cos  —  cos 
—A  cos  <^o 

V?7 


(A.2) 


(A.3) 


(A.4) 


h  =  'J^exp{-jX)(A3  7  J  /J3(f)exp(-A's^)iis)  (A.5) 

V  -4  -4  / 


(A.6) 
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and 


m)  = 


cos  ^  sin 


N 


JV 


,  ^ 
AT 


The  asymptotic  evaluation  of  (A.5)  yields 

I3  ~  — 2j7riVsgn(a)cos(^oexp(— jJ*f)F(|a|\/X) 
-i-exp(-jX)i 


/2j7r 

f  N  cos  ^0 

sin^  \ 

/  X 

\  \/2^ 

cos^  -  cos^y 

(A.7) 


(A.8) 
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Appendix  B 

Evaluation  of  Integrals 


In  this  appendix  the  integrals 

’exi>{—jku  sin^  /3o(l  —  A«)) 


OO 

F,=j 


y/u 


du 


(B.l) 


and 


^2  =  j  F{y/2kusin^o\cos^\)exi>(—jkus\n^0o{l  —  n))du 


(B.2) 


are  calculated.  It  is  noted  that  n<\.  Thus,  applying  the  substitution  F  =  tifcsin^/!9o(l  — 
and  assuming  that  sin^  /?o(l  ~  7^  O5  the  integral  Fi  becomes 


Fi  = 


^A:(l-A.)sin^o^^ 


OO 

I  exp(— di 


(B.3) 


which  is  expressed  in  terms  of  the  modified  Fresnel  function  F  (3,22) 
y/jk{l  -  fi)  sin/9o 

where  L  =  kl^  sin^  For  sin^  ;0o(l  —  A*)  =  0  the  integral  Fi  is  singular. 
The  integral  F2  can  be  rewritten  as 


(B.4) 


F, 


rr  OO  OO 

=  J—  J  exip{jku  sin^  ^oif^  +  cos  <f>o))  j  exp{—jF)  dt  du 

y/2ku  sin  /3o  |  cos  ^  | 


(B.5) 
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which  for  sin^  /?o(^  +  cos  ^o)  ^  0  is  evaluated  using  integration  by  parts 


exp(jku  sin^  +  cos  ^o)) 

j  k  sin^  ^o{n  +  cos  ^o) 


00 

j  exp(-jt^)dt 


fi 


y/ikusin  Po  |  cos  ^  | 


\/fcsin/?o|cos  ^ 


\/2« 

—  exp{jL{fi  +  cos  <^o)) 
■s/jlrk  sin^  l3o{fi  +  cos  ^o) 

Icos^l 


^  '  du 


OO 

J  exj>(—jt^)  dt 


\/2L|cos^| 

F^. 


y/2j'Kk  sin  ^q{ii  +  cos  ^o) 

Using  the  expression  for  Fi  in  (B.4),  the  result 

+  cos^o)  0. 


F,= 


jksia^  +  cos 
is  obtained.  This  expression  has  a  finite  limit  for  sin^ 


(B.6) 


(B.7) 
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MISSION 

OF 

ROME  LABORATORY 


Mission.  The  mission  of  Rome  Laboratory  is  to  advance  the  science  and 
technologies  of  command,  control,  communications  and  intelligence  and  to 
transition  them  into  systems  to  meet  customer  needs.  To  achieve  this, 
Rome  Lab: 


a.  Conducts  vigorous  research,  development  and  test  programs  in  all 
applicable  technologies; 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportability; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Materiel 
Command  product  centers  and  other  Air  Force  organizations: 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence,  reliability 
science,  electro-magnetic  technology,  photonics,  signal  processing,  and 
computational  science. 

The  thrust  areas  of  technical  competence  include:  Surveillance, 
Communications,  Command  and  Control,  Intelligence,  Signal  Processing, 
Computer  Science  and  Technology.  Electromagnetic  Technology, 
Photonics  and  Reliability  Sciences. 


